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ABSTRACT 

We determine the velocity correlation function, pairwise peculiar velocity 
difference, and root-mean-square pairwise peculiar velocity dispersion of rich 
clusters of galaxies, as a function of pair separation, for three cosmological mod- 
els: = 1 and = 0.3 CDM, and O = 0.3 PBI models (all flat and COBE- 
normalized). We find that close cluster pairs, with separation r < 10/i _1 Mpc, 
exhibit strong attractive peculiar velocities in all models; the cluster pairwise 
velocities depend sensitively on the model. The mean pairwise attractive veloc- 
ity of clusters on 5/i _1 Mpc scale ranges from ~ 1700km s _1 for = 1 CDM, 
to ~ 1000km s" 1 for PBI, to ~ 700km s" 1 for O = 0.3 CDM. The small-scale 
pairwise velocities depend also on cluster mass: richer, more massive clusters 
exhibit stronger attractive velocities than less massive clusters. On large scales, 
from ~ 20 to 200/i _1 Mpc, the cluster peculiar velocities are increasingly domi- 
nated by bulk and random motions; they are independent of cluster mass. The 
cluster velocity correlation function, which reflects the bulk motion minus the 
relative motion of pairs, is negative on small scales for = 1 and O = 0.3 
CDM, and positive for PBI; this indicates stronger pairwise motion than bulk 
motion on small scales for CDM, and relatively larger bulk motions for PBI. 
The cluster velocity correlation function is positive on very large scales, from 
r ~ 10/i _1 Mpc to ~ 200/i _1 Mpc, for all models. These positive correlations, 
which decrease monotonically with scale, indicate significant bulk motions of 
clusters up to ~ 200/i _1 Mpc. The strong dependence of the cluster velocity 
functions on models, especially at small separations, makes them useful tools in 
constraining cosmological models when compared with observations. 

Cosmology: large-scale structure of Universe - Cosmology: theory - galaxies: 
clustering - galaxies: clusters of 

1. INTRODUCTION 

Observations of large-scale motions in the universe indicate that the Hubble 
expansion is considerably perturbed on scales up to at least 60/i _1 Mpc (Rubin 
et al. 1976; Collins et al. 1986; Dressier et al. 1987; Burstein et al. 1987; Suther- 
land 1988; Lynden-Bell et al. 1988; Faber et al. 1989; Mould et al. 1991; Strauss 
et a/.1992; Mathewson et al. 1992). The perturbations are believed to be caused 
by the gravitational growth of cosmic structures. The observed peculiar motions 
can therefore place critical constraints on cosmological models. One of the most 
important quantitative statistical tools used in contrasting models with observa- 
tions is the peculiar velocity correlation function: comparing the strength and 
the scale on which galaxy pair motions are correlated with each other (Peebles 
1987, 1993; Gorski et a/.1989; Groth et al. 1989). Almost all models and obser- 
vations so far have concentrated on studying the peculiar velocities of galaxies. 
However, uncertainties in the observational determination of galaxy peculiar ve- 
locities, which are proportional to the galaxy distances, do not allow accurate 
velocity determinations except for the nearest galaxies. 

The velocities of rich clusters of galaxies can be determined more accurately 
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than the velocities of individual galaxies since the distance to each cluster can 
be obtained by averaging the measured distances of a large number of member 
galaxies. The Brightest Cluster Galaxy (BCG) distance method also appears to 
be promising (Lauer & Postman 1994). In addition, cluster peculiar velocities 
can, in principle, be determined to large distances using another, independent 
method of the Sunyaev-Zel'dovich (1980) effect: measuring the Doppler shift of 
the scattered cosmic microwave radiation (CMR) caused by the peculiar motion 
of the cluster relative to the CMR frame. While the effect is small, current 
observations are beginning to use this method to measure cluster motions. Clus- 
ter motions could therefore provide a major (and economical) tool in probing 
the large-scale peculiar velocity field (Bahcall, Gramann & Cen 1994; Strauss et 
a/.1994). 

In this letter we determine three quantitative statistics of the peculiar mo- 
tions of rich clusters of galaxies in different cosmologies: 1) the velocity cor- 
relation function of clusters; 2) the pairwise velocity difference of clusters as a 
function of cluster pair separation; and 3) the root-mean-square pairwise velocity 
dispersion of clusters as a function of pair separation. We also investigate the 
dependence of the above functions on cluster richness. 

2. MODEL SIMULATIONS 

A Particle-Mesh code (Cen 1992) with a box size of 800^ -1 Mpc is used to 
simulate the evolution of dark matter in different cosmological models. A large 
simulation box is needed in order to 1) ensure that contributions to velocities 
from waves larger than the box size are small, and 2) minimize uncertainties due 
to fluctuations in the small number of large waves. Each simulation has 500 3 
cells and 250 3 = 10 7 ' 2 particles. The nominal spatial resolution is 1.6/i _1 Mpc. 

Three cosmological model simulations are performed: Standard Cold Dark- 
Matter (SCDM) model (O = 1, h = 0.5, a 8 = 1.05), low-density CDM (LCDM) 
model (O = 0.3, A = 0.7, h = 0.67, erg = 0.67), and Primeval Baryonic Isocur- 
vature (PBI) model (O = 0.3, A = 0.7, h = 0.5, <j 8 = 1.02). The models are 
normalized to the CMR anisotropy measured by COBE (Smoot et al. 1992). 

Clusters are selected in each simulation using an adaptive, friends-of-friends 
linkage algorithm described by Suto, Cen & Ostriker (1992) and Bahcall & Cen 
(1992). The cluster mass threshold, within r = 1.5/i _1 Mpc of the cluster cen- 
ter, is selected to correspond to a number density of clusters comparable to the 
observed density of rich (R > 1) clusters, n c i ~ 6 X 10 -6 /i 3 Mpc -3 (Bahcall & 
Cen 1993). A total of ~ 3000 rich clusters of galaxies are obtained in each of the 
800/i _1 Mpc simulation models. The three-dimensional peculiar velocity of each 
cluster, relative to the co-moving cosmic background frame, is determined from 
the simulation; these velocities are used in the analyses described below. In order 
to investigate the dependence of the results on cluster richness (i.e., mass), we 
carry out the velocity analysis independently for rich (R > 1) clusters, for the less 
massive R > clusters with a mean density of n cl (R > 0) « 1.3 x 10" 5 /i 3 Mpc- 3 , 
and for lower threshold groups with n gr % 10 -4 /i 3 Mpc -3 (Bahcal & Cen 1993); 



3 



these comparisons are discussed in sections 3 and 4. 

The sensitivity of the results to the resolution of the simulation is tested 
by comparing simulations with 0.8 and 1.6/i _1 Mpc nominal resolutions (both 
using the same 400/i _1 Mpc box). The resulting cluster velocity distributions are 
consistent to better than 60km s _1 . The simulation results are also compared 
with expectations from linear theory (Gramann et al. 1994). 

3. VELOCITY CORRELATION FUNCTION OF CLUSTERS 

The velocity correlation function of rich clusters, ^ v (r) = sign| < vi(fi) ■ 
V2{ri + r) > I 1 / 2 = sign|| < \vi(fi) + v 2 {ri + r)] 2 — \vi(fi) — v 2 {ri + r)] 2 > j 1 / 2 
(Peebles 1987, 1993; Groth et al. 1989), is determined as a function of cluster 
pair separation, r, for all models. \& v {r) measures the difference of two terms — 
the bulk motion less the relative motion of cluster pairs. ^ v is positive if clusters 
in the pairs move in the same direction, and is negative if they move in opposite 
directions. A random distribution of cluster motions yields ^ v ~ 0; bulk-motion 
exhibits ^ v > 0, and attracting or repulsing cluster pairs have ty v < 0. 

Figure 1 shows ^ v (r) for rich (R > 1) clusters of galaxies. It reveals negative 
correlations for the two CDM models and positive correlations for PBI on small 
scales (r < 10/i _1 Mpc), and positive correlations on larger scales (from r ~ 
10/i _1 Mpc to ~ 200/i _1 Mpc) for all models. The small scale negative correlations 
are strongest in SCDM, with ^> v (bh~ 1 Mpc) ~ —700km s _1 . On large scales, 
r > 20/i _1 Mpc, the positive correlation function decreases monotonically with 
scale. Table 1 summarizes the rich cluster velocity correlations. 

What are the implications of the above velocity correlation function? The 
negative correlations on small scales (r < 10/i _1 Mpc) imply that the relative 
motion of close rich cluster pairs dominates over their bulk motion. As we show 
in §4.1, close cluster pairs mostly approach each other, as expected. The cluster 
attractive velocities are largest in SCDM. The PBI model yields positive correla- 
tions even at these small separations, suggesting that, in addition to the cluster 
attractive motions (§4.1), large bulk motions (due to the large-scale bump in 
the power spectrum of PBI) contribute significantly (positively) to ty v > 0. On 
larger scales, r > 10/i _1 Mpc, the cluster velocity correlations become positive, 
and remain positive to very large scales (r > 200/i _1 Mpc). This reflects the 
transition from the local gravitational attraction regime for close cluster pairs 
to the large scale regime where the velocity correlations are dominated by bulk 
motions. The PBI model exhibits the strongest positive velocity correlations on 
scales up to r ~ 10 2 /i _1 Mpc, where it drops rapidly due to its steep slope of the 
power spectrum on large scales (r > 400/i _1 Mpc). The positive correlations de- 
crease monotonically with scale for all models for r > 20/i _1 Mpc as expected for 
large-scale bulk motions. The results are consistent with the detailed large-scale 
cluster motions studies by Gramann et al. (1994) for these models. Comparison 
with expectation from linear theory on large scales shows consistency with the 
simulations of the low-density models, and a slight overestimate of the simulated 
bulk-flow of SCDM by ~ 80km s _1 (Gramann et al. 1994). The positive cluster 
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Figure 1. Velocity correlation function of rich (R ^ 1) clusters of galaxies for the three models 
(§3). 



velocity correlations on large scales suggest a similar behavior for the galaxies. 

The velocity correlation function on small scales depends on cluster mass, or 
richness. We find that poorer, less massive clusters exhibit somewhat weaker at- 
tractive correlations on small scales than massive clusters; the large scale velocity 
correlations, however, are independent of cluster mass. The richness dependence 
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Figure 2. Dependence of the cluster velocity correlation function on cluster richness (R > 1, 
R > clusters and groups), for = 1 (top) and = 0.3 (bottom) CDM. (The large-scale 
= 1 CDM velocities are overestimated in comparison with linear theory by ~ 80km S 
(§3). 



of the cluster velocity correlation function is presented in Figure 2, where we 
compare the velocity correlation function of the rich R > 1 clusters, with that 
of poorer R > clusters, and of groups of galaxies (see §2). The results are 
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presented for SCDM and LCDM; PBI yields trends similar to those of LCDM. 
Since the small scale correlations are mostly due to attractive forces between 
close cluster pairs, we find that massive clusters attract each other more strongly 
than less massive clusters. The same trend is observed in the cluster spatial cor- 
relation function where more massive rich clusters are correlated more strongly 
than poorer clusters (Bahcall & Soneira 1983; Bahcall & West 1992). 

4. PAIRWISE VELOCITIES OF CLUSTERS 
4-1 Velocity Difference of Cluster Pairs 

The mean peculiar velocity difference of cluster pairs along the radius vector 
r = r*2 — r*i joining the pair, Vi2{r) = — ((^2 — £1) • r*/r), is shown in Figure 3; 
the brackets indicate an average over all pairs at the same separation. The 
pairwise peculiar velocity difference Vi2{r) is positive for attracting pairs, it is 
negative for repelling pairs, and it is zero for a random distribution of peculiar 
velocities. Bulk-flow motions do not contribute to v\2- The results show large 
positive pairwise velocities at small separations, v\2(r < 10/i _1 Mpc) ~ 500 to 
1700km s _1 , which decrease monotonically to zero on large scales. All models 
exhibit the same behavior but with different amplitudes (Figure 3). 

The strong mean pairwise cluster velocities for r < 10/i _1 Mpc implies a 
strong mean attraction of cluster pairs on these scales. The attraction is strongest 
for SCDM, with v 12 (r ~ 5/i _1 Mpc) ~ 1700km s -1 , and weakest for LCDM, with 
^12 ~ 5/i _1 Mpc) ~ 700km s _1 . The PBI model yields intermediate values. 
The results are summarized in Table 1. The strong attractive motions in SCDM 
are partly due to the larger cluster masses in this model (larger than observed, 
and larger than LCDM and PBI; Bahcall & Cen 1992). 

The monotonic decrease of ^i2(r) with r reflects the relative decrease of the 
cluster gravitational attraction and the increased contribution of random and 
bulk velocities at large scales. The random and bulk velocities finally dominate 
on large scales, r > 50 — 100/i _1 Mpc, where V12 ~ 0. The amplitude of the cluster 
pairwise velocity function on small scales depends on cluster richness, or mass, 
similar to the richness dependence of the cluster velocity correlation function 
(§3); more massive clusters attract each other more strongly than less massive 
clusters on small scales (r < 10/i _1 Mpc). 

4-2 RMS Pairwise Velocity Dispersion of Clusters 

The ID rms pairwise peculiar velocity dispersion of clusters, a\® = ((^1 — 
^2) 2 } 1//2 /"\/3, is presented as a function of pair separation, r, for all models, in 
Figure 4. Unlike Vi2{r) discussed in §4.1, o\® does not determine the mean 
relative direction of the pairs motion; rather, it provides the average magnitude 
of the rms velocity dispersion of cluster pairs. 

The cluster velocity dispersion exhibits large pairwise velocities on small 
scales (~ 400 to 1300km s _1 for r < 10/i _1 Mpc, depending on the model); 
the velocity dispersion decreases nearly monotonically with scale, similar to the 
u i2(0 function, and reaches a large-scale velocity dispersion of ~ 300km s _1 
(LCDM) and ~ 600km s" 1 (SCDM and PBI) on ~ 100fc -1 Mpc scale. All models 
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Figure 3. Cluster pairwise velocity difference as a function of pair separation for the three 
mode Is (§4.1). (R > I clusters). 

exhibit a similar behavior, but have different amplitudes (Figure 3 and Table 1). 
The broad minimum near 30/i _1 Mpc reflects the transition beyond which the 
bulk motions become small relative to the overall velocities. 

The large velocity dispersion on small scales reflects the strong mean attrac- 
tion of clusters on these scales, as implied from V12 (§4.1). The ratio v\2jo\^ 
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Figure 4. Cluster ID pairwise rms velocity dispersion as a function of pair separation for the 
three mode Is (§4.2). (R > 1 clusters). 

is nearly \/3 at r ~ 5/i _1 Mpc, as expected if the cluster pairwise motions on 
these scales are mostly one-dimensional attractive motions along the line joining 
the cluster pair; in that case V12 ~ cf^ ~ y/S&l^ , as observed. This attractive 
motion is strongest for SCDM, with cr^(5/i _1 Mpc) ~ 1000km s _1 , and is lower 
for PBI (~ 650km s _1 ) and LCDM (~ 500km s _1 ). The results are summarized 
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in Table 1. The large difference in the cluster pairwise velocity dispersion on 
these scales makes this parameter an important tool in constraining cosmologi- 
cal models when compared with observations. 

The decrease of the cluster velocity dispersion with scale to r ~ 
and its flattening at a constant positive dispersion on large scales, reflects the 
transition from the gravitational attraction regime of close cluster pairs, to the 
overall rms pairwise velocities of clusters on large scales (Fig. 4, Table 1). 
Cluster pair velocities are mostly uncorrelated on large scales. In this case 
o\® ~ V^u^ s , where vj,® s is the ID rms velocity of clusters. This relation 
is indeed satisfied: v].^ s of rich clusters is 446, 254, and 411km s _1 for SCDM, 
LCDM, and PBI models, respectively (Gramann et al. 1994); these values are 
~ v^2 lower than the large-scale a\® (100/i _1 Mpc) dispersions listed in Table 1. 

The ratio v\2jo\® decreases monotonically with scale from \/Z on small 
scales to zero at r ~ 100/i _1 Mpc. The trend is similar in all models. This 
behavior is expected if the cluster pairwise motions change from one-dimensional 
attraction on small scales to bulk and finally random motions on large scales 
[where v\2 decreases to zero, and a\® to a constant, as observed]. 

The amplitude of the cluster pairwise velocity dispersion on small scales (r < 
10/i _1 Mpc) depends on cluster mass; more massive clusters yield larger pairwise 
velocity dispersions, especially for SCDM. The richness (mass) dependence is 
consistent with that observed for the velocity correlation function (§3) and Vi2{r) 
(§4.1); it suggests stronger attractive motions for more massive clusters on small 
scales. The large scale motions do not depend on cluster mass. 

5. CONCLUSIONS 

We determine the velocity correlation function of rich clusters of galaxies, 
and the pairwise peculiar velocity difference and rms pairwise peculiar velocity 
dispersion of clusters for three cosmological models: = 1 and = 0.3 CDM, 
and O = 0.3 PBI. Our principal conclusions are listed below. 

1. All models show that close cluster pairs (r < 10/i _1 Mpc) exhibit strong 
attractive velocities. The mean pairwise attractive velocity on 5/i _1 Mpc scale 
ranges from ui2(5/i _1 Mpc) ~ 1700km s _1 for = 1 CDM to ~ 700km s _1 
for O = 0.3 CDM. The pairwise cluster velocity function decreases monoton- 
ically with scale, reflecting the relative decrease of the pairwise cluster grav- 
itational attraction and the relative increase of random and bulk motions on 
large scales. The random and bulk motions finally dominate on large scales, 
r ~ 50 — 100/i _1 Mpc, where v\2(r) ~ 0. Massive clusters exhibit somewhat 
stronger attractive velocities than less massive clusters on small scales; on large 
scales, where the cluster velocities are dominated by bulk and random motions, 
the velocities are independent of cluster richness. 

2. The rich cluster velocity correlation function reveals positive correlations 
to very large scales, from r ~ 10/i _1 Mpc to 200/i _1 Mpc, for all models (Fig 1). 
The positive correlations decrease monotonically with scale for r > 20/i _1 Mpc, 
approaching a nearly random distribution of cluster velocities at r > 200/i _1 Mpc. 
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The positive correlations on ~ 20 — 200/i _1 Mpc scale reflect bulk motions of 
clusters, produced by large-scale gravitational perturbations. On small scales, 
r < 10/i _1 Mpc, the velocity correlations are negative for = 1 and O = 0.3 
CDM, reflecting strong mean attractive motions for close cluster pairs, especially 
for = 1. On scales up to 10 2 /i _1 Mpc, PBI exhibits the largest positive velocity 
correlations of clusters, reflecting strong bulk motions. 

3. The ratio v^jo^ decreases with scale from ~ \p3> on small scales (~ 
5/i _1 Mpc) to zero on large scales. This demonstrates that clusters in close pairs 
move primarily along the line joining the pair while pairs with large separations 
move mostly randomly. 

4. The rich cluster pairwise velocity functions and velocity correlation func- 
tion discussed above depend strongly on the cosmological model (shape and am- 
plitude of the power spectrum and O). These functions (or some variants of 
them) can therefore serve as a potentially powerful tool in constraining cosmo- 
logical models when compared with observations. 

It is a pleasure to acknowledge NCSA for use Convex-3880 supercomputer. 
Support by NASA grant NAGW-2448, NSF grants AST90-20506, AST91-08103, 
AST93-15368 and HPCC ASC-9318185 are gratefully acknowledged. 

REFERENCES 

Bahcall, N.A., & Cen, R. 1992, ApJ(Lett), 398, L81 

Bahcall, N.A., & Cen, R. 1993, ApJ(Lett), 407, L49 

Bahcall, N.A., & Soneira, R.M. 1983, ApJ, 270, 20 

Bahcall, N.A., & West, M.J. 1992, ApJ, 392, 419 

Bahcall, N.A., Gramann, M., & Cen, R. 1994, ApJ, Nov 20 

Burstein, D., Davies, R.L., Dressier, A., Faber, S.M., Lynden-Bell, D., Terlerich, 
R., & Wegner, G. 1987, in "Galaxy Distances and Deviations from Universal 
Expansions", ed. B. Madore and B. Tully (Dordrecht:Reild), pl23 

Cen, R. 1992, ApJS, 78, 341 

Collins, A., Joseph, R.D., and Robertson, N.A. 1986, Nature, 320, 506 

Dressier, A., Faber, S.M., Burstein, D., Davis, R.L., Lynden-Bell, D., Terlevich, 
R.J., and Wegner, G. 1987, ApJ(Letters), 313, L37 

Faber, S.M., Wegner, G., Burstein, D., Davies, R.L., Dressier, A., Lynden-Bell, 
D., Terlerich, R.J. 1989, ApJS, 69, 763 

Gorski, K.M., Davis, M., Strauss, M.A., White, S.D.M., and Yahil, A. 1989, ApJ, 
344, 1 

Gramann, M., Bahcall, N.A, Cen, R., & Gott, J.R. 1994, ApJ, submitted 
Groth, E.J., Juszkiewicz, R., & Ostriker, J. P., 1989, ApJ, 346, 558 



- 11 - 



Lauer, T.R., & Postman, M. 1994, ApJ, 425, 418 

Lynden-Bell, D., Faber, S.M., Burstein, D., Davis, R.L., Dressier, A., Terlevich, 
R.J., and Wegner, G. 1988, ApJ, 326, 19 

Mathewson, D.S., Ford, V.L., & Buchhorn, M. 1992, ApJS, 81, 413 

Mould, J., Staveley-Smith, L., Schommer, R.A., Bothum, G., Hall, P.J., Han, 
M., Huchra, J. P., Roth, J., Walsh, W., & Wright, A.E. 1991, ApJ, 383, 467 

Peebles, P.J.E. 1987, Nature, 327, 210 

Peebles, P.J.E. 1993, "Principles of Physical Cosmology", (Princeton University 
Press, Princeton, NJ) 

Rubin, V.C., Thonnard, N., Ford, W.K., & Roberts, M.S. 1976, AJ, 81, 719 

Smoot, G.F., et a/.1992, ApJ(Letters), 396, LI 

Strauss, M., Davis, M., Yahil, A., & Huchra, J. P. 1992, ApJ, 385, 421 

Strauss, M., Cen, R., Ostriker, J. P., Lauer, T., & Postman, M. 1994, ApJ, 
submitted 

Sunyaev, R., & Zel'dovich, Ya.B. 1980, MNRAS, 190, 413 

Sutherland, W. 1988, MNRAS, 234, 159 

Suto, Y., Cen, R., & Ostriker, J. P. 1992, ApJ, 395, 1 

TABLE 1 

VELOCITY CORRELATIONS AND PAIRWISE VELOCITIES 
OF RICH (R > 1) CLUSTERS* 

Model t/>„(5) ip v (20) tp v (100) v 12 (5) <rj^(5) (7^(100) 

CDM O = 1.0 -692 426 235 1698 1070 595 

CDM O = 0.3 -87 275 162 714 487 327 

PBI O = 0.3 451 554 225 962 646 555 



* The velocities are in km s _1 . The scale at which each parameter is listed, indicated 
inside the parentheses, is in units of /i _1 Mpc. The statistical pair uncertainties are 
approximately 15% at 5/i _1 Mpc, 5% at 20/i _1 Mpc, and lower at large scales. 



